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Abstract. Photoproduction of neutral and charged pions off nucleons and deuterium 
has been precisely calculated in baryon chiral perturbation theory. I review the predic- 
tions in light of the accurate data that have become available over the last few years. 
Some progress in the description of neutral pion electroproduction off protons is also 
discussed. 



1 Introduction 

With the advent of CW machines, pion production by real and virtual photons 
has become a major testing ground for predictions based on nuclcon chiral dy- 
namics. In particular over the last years there has been considerable activity to 
precisely measure pion photo-and electroproduction in the threshold region at 
various laboratories, like e.g. MAMI, SAL and NIKHEF. On the other hand, 
refined calculations within heavy baryon chiral perturbation theory (HBCHPT) 
have been performed. Pion photo-and electroproduction is not only interesting 
per se but also allows, as we will see in the following, to get informations on 
other processes like for example irN scattering. In this respect the use of CHPT, 
see for example Gasser et al. (1984), is particularly advantageous since it is a 
method for solving QCD at low energy which links different processes in a model- 
independent fashion thus allowing for a deeper understanding of the underlying 
dynamics. Another very important problem which can be addressed in dealing 
with these processes is isospin symmetry. It has always been one of the major 
goals in nuclear physics to understand isospin symmetry violation related to the 
light quark mass difference m u — and virtual photon effects. Although the 
light quark mass ratio deviates strongly from unity, md/m u ~ 2, see Gasser et 
al. (1975) and one could expect sizeable isospin violation, such effects are effec- 
tively masked since — m u <C A, with A the scale of the strong interactions 
(which can be chosen to be AttF^ or 1 GeV or the mass of the p). To assess 
the isospin violation through quark mass differences, precise measurements and 
accurate calculations are mandatory. As pointed out by Weinberg (1977) long 
time ago, systems involving nucleons can exhibit such effects to leading order 
in contrast to the suppression in purely pionic processes due to G-parity. Thus 
pion photoproduction is well suited for such investigations since a large body of 
(precise) data exists for various isospin channels. 

Here I will report on the progress made since the MIT workshop, see MIT 
proceedings (1994), on the calculation of charged and neutral pion photo-and 
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electroproduction near threshold in the framework of HBCHPT. In this frame- 
work the nucleons are treated as very heavy static sources. HBCHPT is a triple 
expansion in external momenta, quark masses and inverse powers of the nucleon 
mass (collectively denoted by the small parameter q). Calculations which will 
be discussed here are 0(q n ) with n=4 for the S-wave and n=3 for the P-waves 
and only take into account isospin breaking effects which are believed to be 
dominant, namely those through the pion mass difference in the loops. 

2 Formal aspects 

2.1 Effective field theory 

The effective Lagrangian which will be needed here consists of the following 
pieces: 

r _ r(2) , r (i) , /.(l) , A 2 ) A3) AA) , AO) A 1 ) (1 \ 

where the index (i) gives the chiral dimension (number of derivative and/or 
meson mass insertions). L^l + is the non- linear a model Lagrangian cou- 
pled to nucleons. The terms from L^ N +L^ N contributing to the single-nucleon 

photoproduction amplitudes are given in Bernard et al. (1996(1)). C^n nas been 
used in the nuclear force calculation of van Kolck (1994) and enters the calcu- 
lation of (neutral) pion photoproduction on deuterium by Beane et al. (1997), 
which will be discussed below. In HBCHPT the nucleon mass term is replaced by 
a string of 1/m suppressed interactions so that n > 2, contains 1/m terms 
as well as counterterms whose coefficients are the famous low-energy constants 
(LECs) . Some of these counterterms cancel the divergences of certain loop dia- 
grams and thus are scale dependent, Cj = C^/i) + ZiL(/i) where /x is the scale 
of dimensional regularization (naturally the physical amplitudes are scale inde- 
pendent). In the following these LECs which are not fixed by chiral symmetry, 
will either be fitted to experimental data or will be obtained using the principle 
of resonance saturation. In the meson sector it can indeed be shown that the 
numerical values of the renormalized LECs L^(fj, = M p ) can be understood to a 
high degree of accuracy from resonance saturation, i.e. they can be expressed in 
terms of resonance masses and coupling constants of the low- lying vector ( V) , 
axial vector (A), scalar (S) and pseudoscalar (P) multiplets, see Ecker et al. 
(1989). In the nucleon sector there exists no proof of this principle, however 

(2) 

it seems to work rather well in the case of C^^ as has been demonstrated in 
Bernard et al. (1997(1)) and can be seen in table 1 which gives the values of 
the seven finite LECs at that order. The first four have been determined (sec- 
ond column) from a best fit to a set of nine subthreshold and threshold nN 
observables that to one-loop order q 3 are given entirely in terms of tree graphs 
including insertions from these LECs and finite loop contributions, but with 

(3) 

none from the 24 new LECs of C\L. The other three can be determined from 
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the strong neutron-proton mass difference (05, which is only relevant in the case 
m u 7^ ma) and from the anomalous magnetic moments of the proton and the 
neutron (c 6 , C7). Note that the scalar mass to coupling constant ratio Ms/y/gs 
needed to saturate the LEC c\ is in perfect agreement with typical ratios ob- 
tained in boson-exchange models of the NN force, where the a-meson models 
the strong pionic correlations coupled to nuclcons. Note also that the impor- 
tant A degree of freedom is included in the determination of the LECs through 
resonance exchange. I will come back to this point in the following. 
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c Ros ran g eg 


Res 


1 


-0.93 ±0.10 


-0.9* 




S 


2 


3.34 ±0.20 


3.9 


2. ..4 


A,R 


3 


-5.29 ±0.25 


-5.3 


-4.5... -5.3 


A,R,S 


4 


3.63 ±0.10 


3.7 


3.1. ..3.7 


A,R,p 


■5 


-0.09 ± 0.01 








6 


5.83 


6.1 




P 


7 


-2.98 


-3.0 




P,LO 



Table 1: Values of the LECs a in GeV -1 for i = 1, ... ,5. The LECs cqj are 
dimensionless. Also given are the central values (cv) and the ranges for the Ci 
from resonance exchange. The * denotes an input quantity. R and S denote the 
Roper and the scalar resonances respectively 



2.2 Threshold pion photo-and electroproduction 

Consider the process j(k)+N(pi) — > ir a (q)+N(p2), with N denoting the nucleon 
(proton or neutron), 7 a real (k 2 = 0) or a virtual (fc 2 < 0) photon and ir a a 
pion of isospin a. The polarization vector of the photon is denoted by e„. In 
the threshold region, the three- momentum q of the pion is small and vanishes 
at threshold. It is therefore advantageous to perform a multipole decomposition 
since at threshold only the S-waves survive and close to threshold one can confine 
oneself to S- and P-waves. The corresponding multipoles are called (E, M, L)i±, 
where E, M, L stands for electric, magnetic and longitudinal (this last one of 
course only comes into play for virtual photons), I = 0, 1, 2, . . . the pion orbital 
angular momentum and the ± refers to the total angular momentum of the 
pion-nucleon system, j = I ± 1/2. These multipoles parametrize the structure 
of the nucleon as probed with low energy photons. Consequently, the T-matrix 
depends on seven (only four survive for real photons) complex multipoles and 
takes the following form: 

T ■ e = ier ■ e (e 0+ + q ■ k P^j + i<r ■ k e • q P 2 + {q x k) ■ e P 3 
+i<r-ke-k ( L + -E 0+ +q-k (P 4 -P 5 -Pt- P 2 ) j +icr ■ qe ■ kP 5 (2) 



4 Veroniquo Bernard 



The quantities P1-5 represent the following combinations of the five P-waves, 
E 1+ , Mi+, Mi_, 

Pi = 3Pi+ + Mi+ - Mi_ , P 2 = 3Pi+ - Mi+ + Mi_ , P 3 = 2Mi+ + Mi_ 

P 4 - 4Li+ + Li_ , P 5 = Li_ - 2L 1+ . (3) 

All these amplitudes are easily calculable within CHPT. They have the conven- 
tional isospin decomposition (to first order in the electromagnetic coupling), 

A(a, u) - #) (a, U )(5 a , 3 + (s, u)± [r a , r 3 ] + A< ' (s, u)r 3 (4) 

There are however, four physical channels, two charged reactions and two neutral 
ones which will be discussed next. Thus there exists a triangle relation relating 
one of the physical amplitudes to the others. This relation should of course only 
hold if isospin is an exact symmetry. It is then clear that it is very important 
to have a very precise determination of the four physical channels in order to 
measure isospin violation. 



3 Charged Pion Photoproduction 

Charged pion photoproduction is a particularly interesting process since it allows 
to 

• investigate the violation of isospin symmetry beyond leading order in elec- 
tromagnetism. For that one of course needs to know extremely precisely the four 
physical photoproduction reactions as discussed previously. 

• determine irN scattering lengths. 

• give a stringent constraint on the much discussed value of the pion-nucleon 
coupling constant via the Goldbcrger-Miyazawa-Oehme (1955) sum rule 
combined with the Panofsky ratio. 

It is well described by the Kroll-Ruderman term which is non-vanishing in the 
chiral limit, 



with fj, = M w + /to and using gl N /4ir = 14.28, e 2 /47r = 1/137.036, to = 928.27 MeV 
and Af 7r+ = 139.57 MeV. In the limit M n = 0, this simplifies to 

E^(ir + n) = -E^(ir-p)=34-W- 3 /M 7r . (6) 

By comparing the numbers in Eq.(5) and Eq.(6) one notices that the kinematical 
corrections which are suppressed by powers of the small parameter /j, ~ 1/7 
are quite substantial for Pg^ r (7r + n). However, there are other corrections which 
are related to pion loop diagrams and higher dimension operators. These have 
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been dealt with in a systematic fashion using heavy baryon CHPT up-to-and- 
including order 0(/z 3 ) in Bernard et al. (1996(2)). In this framework, one has to 
consider pion loop diagrams and local contact terms whose coefficients are the 
LECs. These were estimated by resonance exchange since not enough precise 
data exist to pin them all down. Frozen kaon loops contributes to E^) and 
while p-meson exchange contributes only to E^) and the A, the axial resonance 
and the C!_ 3 (see previous section) to -Eq + \ where E 0+ (ir + n) = \/2(E^ + Eq^ ) 
and £o+(7r-p) = \/2(4+ } - E o+ y )- This 01 

course leads to some uncertainty 
into the result since the resonance parameters are only known within certain 
ranges. Another source of uncertainty comes from the regularization scale. Indeed 
within resonance saturation a spurious mild scale-dependence remains. In the 
calculation A runs in the interval M p < A < m^. To C(^ 3 ) one gets with 
g„ N = 13.4: 

E$ = (-1.6 ± 0.1) • KT 3 /^ , E$ = (21.5 ± 0.4) • 10~ 3 /M„ (7) 

The value of g„N will be discussed in the next section. The chiral expansion is 
rapidly converging in contrast to neutral pion photoproduction: 

E$ = (0 - 1.79 + 0.38 - 0.21) • 10- 3 /M w , 

E$ = (24.01 - 3.57 + 1.38 - 0.29) • lO^/M^ (8) 

where the various contributions to E$ and E ( Q+ ] of C(M") with n=0,l,2,3 have 
been collected. Translating these results into the physical channels one obtains 
the results shown in the second column of table 2. 





CHPT 


DR 


Experiment 




28.2 ±0.6 
-32.7 ±0.6 


28.0 ±0.2 
-31.7 ±0.2 


27.9 ± 0.5, 28.8 ± 0.7, 27.6 ± 0.3 
-31.4 ± 1.3, -32.2 ± 1.2, -31.5 ± 0.8 



Table 2: Predictions and data for the charged pion electric dipole amplitudes. 

Also given in that table are the results of the dispersion theoretical (DR) analysis 
of Hanstein et al. (1997) and the experimental ones. The first two numbers in 
the last column corresponds to rather old data from Burq (1965), Adamovitch 
(1966) and Goldwasser et al. (1964) while the last ones are taken from the 
recent TRIUMF experiment on the inverse reaction tt~p — > jn (see Kovash et 
al. (1997)) and a preliminary SAL analysis of the reaction — > ir + n. The overall 
agreement is quite good though the CHPT predictions lie on the large side of 
the most recent data. Clearly, we need more precise data together with a better 
knowledge of the pion nucleon coupling constant to draw a final conclusion. 

With the determination of the threshold value of Eo+ for the reaction — > 
n~p one can deduce the value of the difference between the isospin 3/2 and the 
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isospin 1/2 nN scattering length through the following relation based on time 
reversal invariance 

(a 3 -ai) 2 =9P^E 0+ . thl ( 1 n^ir-p) (9) 

where P is the experimentally well determined Panofsky ratio (P = u(t:~p — > 
Tr°n)/a(Tr~p — > -yn) = 1.543 ± 0.008) and q 7 and are the CM momenta of 
the photon and neutral pion at the pir~ threshold, respectively. In the third 
column of table 3, the value of the isovector scattering length — bi = (a^ — 
ai)/3 as obtained from Eq.(9) is given. It is compared with a direct CHPT 
calculation of the scattering process irN — > irN (note that in that case the 
result is independent of g^N to leading order, see Bernard et al. (1995)). The 
two values are consistent within the error bars. These results clearly demonstrate 
the importance of looking at different processes at the same time. One certainly 
needs a better determination of the LECs, a better control on the residual scale 
dependence and isospin breaking effects have to be included before one can reach 
any definite conclusion. For comparison are given the Karlsruhe-Helsinki (KH) 
result, see Hohler (1983) and the value obtained by Sigg et al. (1995) from the 
decay width in pionic hydrogen. 





CHPT 

irN -> nN 


CHPT 

Eo+ 


KH 


decay width 




92 ±4 


87.3 ±2.4 


91.3 ±4 


96 ±7 



Table 3: Isovector irN scattering length. 

4 Neutral Pion Photo-production 

4.1 7p — > tv°p 

Neutral pion photoproduction has been a hot topic ever since the Saclay and 
Mainz groups claimed a sizeable deviation from a so-called low energy theorem 
(LET) for the electric dipole amplitude E 0+ derived in 1970. In the MIT proceed- 
ings (1994), I reported on preliminary results which were then confirmed, from 
an HBCHPT calculation by Bernard et al. (1996(1)). I will briefly summarize 
some by now well known facts. 

• E n+ : To 0(q 3 ) in HBCHPT two graphs contribute, the famous rescat- 
tcring graph and the triangle graph. This last important one brings in some 
non-analytical pieces which invalidate the assumptions used to derive the so- 
called LET of the seventies, this one thus just being a LEG (low energy guess) 
(see Ecker and Meifiner (1995)). To 0(q 4 ) one has to take into account some 
more loop diagrams (1/m corrections to the one just discussed and loops with 

(2) 

one vertex or with a propagator from I^-iv) P ms * w0 counterterms. It was 
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Fig. 1. The real part of the electric dipole amplitude in the threshold region. Left 
panel:7p — > Tr°p. The dashed line is the prediction of the incomplete LET and the data 
are from the multipole analysis of Fuchs et al. (1996) (boxes) and of Bergstrom et al. 
(1997) (crosses). Right panel:7n — + ir°n 

shown that the expansion of E + in powers of fi = M^/m is slowly converg- 
ing (E + = -3.45(1 - 1.26 + 0.55) where the M % , Ml and M% contributions are 
given, to be compared with Eq.(8)) and therefore hard to pin down accurately. 
It was concluded that at threshold this multipole cannot provide the best test 
of chiral dynamics as was first conjectured. 

• P-waves: LETs have been derived for P x 2 . In contrast to E n+ these are 
very fastly converging functions of ^. It is thus particularly important to test 
these directly. As we will see polarization measurements have and will be made 
for that purpose. It was also shown that the P-waves are of the same chiral order 
as Eo + , namely they are proportional to \q\ and not \q\\k\ as usually postulated, 
where q and k are the pion and photon cm momenta respectively. P\ t 2 have 
a very week energy dependence, they show no cusp effect and they have very 
small imaginary parts. To 0(q 3 ) P3 is completely dominated by a counterterm b p 
which in the resonance saturation picture is essentially described by the Z\(1232) 
resonance . 

It turned out that new data from the TAPS collaboration (see Fuchs et al. 
(1996)) were released and showed some discrepancies to the previously considered 
best data of Beck et al. (1990). They were analysed in Bernard et al. (1997(2)) 
using the same HBCHPT formalism as in 1996. The total and differential cross 
sections were fitted and E + was then predicted, the result is shown in Fig.l. 
Note that even so the convergence for this particular observable is slow, a CHPT 
calculation to order p 4 allows to understand its energy dependence in the thresh- 
old region. Indeed nice agreement between theory and experiment is obtained as 
confirmed by table 4 which gives the values of Eq + at the ir°p and 7r + n thresh- 
old. For comparison is also shown the result of the dispersion theoretical analysis 
(DR) of Hanstein et al. (1997). The reaction -yp — > n°p was also remeasured at 
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Saskatoon, see Bergstrom et al. (1996). Between tt°p and ir + n thresholds, the 
SAL data are consistent with the ones of Fuchs et al. (1996). For larger energies, 
however, the new SAL data agree with the older Mainz data, see Beck et al. 
(1990). This does not affect the threshold value of E 0+ but rather leads to a 
larger value of b p . The experimental discrepancy remains to be clarified. It was 
pointed out by R. Beck at this workshop that a measurement of the so called 
F asymmetry in a polarized photon (circular) and polarized target (x-dircction) 
7P — ► n°p experiment would allow to determine RcE n+ with a small statistical 
error. 





CHPT 


DR 


Experiment 


E a+ (TT° P thl) 

£ , 0+ (7T + nthr) 


-1.16 

-0.44 


-1.22 
- -0.4 


-1.31 ± 0.8, -1.32 ± 0.05 ± 0.06 
0.4 



Table 4: Predictions and data for the electric dipole amplitude for the neutral 
pion photoproduction off protons at the ir p and ir + n threshold. 

Since the MIT workshop no real progress was made on Imi?o+. This is still 
an important piece of information missing on the experimental side since it is 
related to the change of KcE 0+ through a dispersion relation as well as to irN 
scattering length through unitarity. A recent multipolc analysis of the Mainz 
■yp — > n°p cross section measurements by Bernstein et al. (1997(1)) seems to 
indicate some sensitivity to the assumed energy dependence of the P-wave mul- 
tipolcs. As already emphasized in the MIT proceedings (1994) a measurement 
of the polarized target asymmetry T (T is sensitive to a linear combination of 
P-wave multipoles times InxEo+) is awaited (see also Bernstein (1997(2)) for a 
discussion of this topic). There are actually two proposals for doing a polarized 
target (y-direction) 7p — > ir°p experiment, one at DUKE (B. Norum) and one 
at MAMI (A. Bernstein). 

The present status of what is known from unpolarized target experiment 
concerning the P-waves is summarized in table 5. There are given the three P- 
waves pi, P23 with p| 3 = {p\ + p\)/2 and e i+ where the small letters refer to 
the reduced multipoles Pi = Pi\k\\q\ where Pj is defined in Eq.(2). I kept here 
the standard assumed behaviour of the P-waves though as pointed out before, 
it is not the proper one. Experimental numbers are from Fuchs et al. (1996) 
(first number for p\ and P23 ) and Bergstrom et al. (1997) (second number for 
pi and P23 and number for ei+). Theory and experiment agree quite nicely for 
pi and p23- Note that there is no prediction for p 2 3 in the case of CHPT since 
as pointed out previously, the value of the p 3 mutipole mainly comes from a 
counterterm fitted to experiment. Still it is important to notice that the value of 
bp is in excellent agreement with the resonance exchange estimate. In the absence 
of polarization data, a unique separation of the three P-wave multipoles is not 
possible. Thus in Bergstrom (1994) milder constraints provided by the Virginia 
Polytechnic Institute multipole analysis, and certain theoretical considerations 
have been applied to effect the separation, thus leading to a value for e\ + . The 
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one given in table 5 comes from a refined analysis of pion angular distributions, 
however the error bar is so large that it is hard to conclude anything at present. 
Also note that the CHPT calculation was done to 0(q 3 ) which means that only 
the first term in the fx expansion of e 1+ is known, ft is certainly necessary to 
go one order higher to see how big the first correction is. Work by the BKM 
collaboration in this direction is in progress. 





CHPT 


DR 


Experiment 


Pi 


10.33 


10.52 


10.02 ±0.15, 10.26 ±0.10 


P23 


fit to exp 


10.85 


11.44 ±0.09, 11.62 ±0.08 


ei 


-0.11 


-0.15 


-0.25 ±0.16 (from unpol.) 



Table 5: Predictions and data for the reduced p- waves multipolcs in units of 
\q\\k\10^/M 3 + 

What is certainly needed is a measurement of the photon asymmetry S in a 
polarized photon experiment. Indeed this quantity defined as 

m . 

e 7 da + da" 

where e 7 is the degree of linear photon polarization, is proportional to P^- 

Thus measuring S allows to test the LET for P 2 and also knowing the multipolcs 
Pi and P23 from unpolarized target to have an unambiguous determination of the 
three multipoles Mi+, Mi_ and E\ + . Data taking was completed in Mainz by 
Ahrens et al. (1994) in the threshold region (E y = 145 - 190 MeV). Preliminary 
results were presented by R. Beck at this workshop. The preliminary measured 
value of £ of (8 ±5%) for 6 = 80° - 100° and E 1 = 155MeV is in good agreement 
with the CHPT prediction of 10%. 

In the calculations presented here the 4(1232) resonance is not treated as an 
explicit degree of freedom. It enters in the counterterms through the principle 
of resonance saturation. It was argued by Jenkins and Manohar (1991) that this 
resonance should be explicitely taken into account since its mass is very close to 
the nucleon mass (A — vtia — m n ~ 3F W ) and the couplings of the N A system 
to pions and photons are very strong, e.g. g n NA ~ ZgNN-x- Recently, Hemmert 
et al. (1997(1)) proposed a systematic way of including the Z\(1232) based on 
an effective Lagrangian of the type £ c e [U, N, A] which has a systematic " small 
scale expansion" in terms of three small parameters (collectively denoted as e) , 
E n /A, M n /A, A/ A, with A e [M p ,m N A^F^}. The method has been applied in 
particular to the calculation of Eq + in neutral pion photoproduction off protons 
at threshold, see Hemmert et al. (1997(2)). The result obtained is quite similar 
to the one of Bernard et al. (1997(2)), meaning that in that case no explicit 



(10) 
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A is needed which is quite reasonable. The situation might be different for the 
P-waves. There one indeed expect larger sensitivity to the A. The inclusion of 
the A as an explicit degree of freedom for these multipolcs is in progress, see 
Bernard et al. (1997(3)). 

In all the calculations done so far the pion coupling constant g^NN was fixed 
to the KH value of 13.4. As pointed out by various persons, the chiral predictions 
thus obtained for E£+ n , E^ n and p are consistently somewhat too big as 
compared to experiment. Decreasing g^NN to 13.06 ± 0.15 one gets very good 
agreement with experiment for this set of observables. This value turns out to be 
consistent with various recent determinations, see the proceedings of the Seventh 
International Symposium on Meson-Nucleon Physics and the Structure of the 
Nucleon, Vancouver (1997). Of course, one has to check what happens for all 
quantities. Furthermore, all isospin breaking effects have to be included before 
one can draw any conclusion. This is one of the major directions which remains 
to be explored on the theoretical side. Work in this direction has started since 
the effective Lagrangian including virtual photons has already been constructed 
by Meifiner and Steininger (1997). 

4.2 -fn — ► 7r°n 

As already discussed, it is very important to have predictions for the reaction 
■yn — > TT°n which is experimentally very difficult to assess. Having fitted the LECs 
to the proton case Bernard et al. (1996(1)) have determined Eq + (it n). A much 
better convergence is obtained due to the fact that the famous triangle diagram 
already appears to first order which is already 0(M%). The first corrections are 
<~ 30%. The result is shown on the right panel of Fig.l. The amplitude clearly 
exhibits the unitary cusp due to the opening of the secondary threshold -/n — > 
n~p — > TT°n. It is sizeably larger in magnitude than for the proton (compare 
with the left panel). This is however not the case at threshold within dispersion 
relations. These indeed tend to give values of the same size, E 0+ (TT°n) = 1.19 
to be compared with the CHPT result E 0+ (ir°n) = 2.13 (note that the DR 
treatment for the neutral channels is less stable than for the charged ones, see 
discussion in Hanstein et al. (1997)). 

5 Neutral pion photoproduction off the deuteron 

The elementary neutron amplitude jn — > n a n can only be inferred indirectly 
from reactions involving few-nucleon systems like the deuteron or 3 He. Beane 
et al. (1997) concentrated on coherent neutral pion production off deuterium in 
the threshold region and studied the sensitivity of the deuteron electric dipole 
amplitude Ed to the elementary neutron amplitude £?o + ™. To tackle this problem 
it is necessary to use a framework which allows one to systematically include 
and order the various contributions arising from single and multiple scattering 
processes. Following the idea of Weinberg (1992) one calculates matrix element 
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of the type < <Pd|/C|^ > by using deuteron wave functions ^ obtained from 
accurate phenomcnological NN potentials and chirally expands the kernel K,. 
Using a large variety of these potentials allows one to assess to which degree of 
accuracy one is sensitive to the chiral symmetry constraints used in determining 
the irreducible scattering kernel. 

The deuteron S-wave multipole is defined in a similar way as the nucleon one: 
T = 2iJ ■ eEd with J = L + S the deuteron total angular momentum. The slope 
of the differential cross section at threshold takes the form: |k|/|q|dcr/dl7| t hr = 
8/3\Ed\ 2 - There are two types of contributions to E^. 

• single scattering contribution 

The single scattering contribution is given by all diagrams where the photon 
and the pion are absorbed and emitted, respectively, from one nucleon with 
the second nucleon acting as a spectator (the so-called impulse approximation), 
leading to 

= IT^ (<P + <n) I d3p 0}(P) 6 ' J Mp - k/2) 

-Afc-| d 3 pp±(pf p + P?' n )cf>}(p)e-JMp-k/2)} , (12) 

evaluated at the threshold value |k| = fc thr = M n o - M 2 /(2m d ) = 130.1 MeV 
and with J = (<ri + <x 2 )/2. A number of remarks concerning Eq.(12) are in 
order. It is important to differentiate between the ir°d and the ir°N (N = p, n) 
center-of-mass (COM) systems. At threshold in the former, the pion is produced 
at rest, it has, however, a small three-momentum in the latter, see Koch et al. 
(1977). Consequently, one has a single-nucleon P-wave contribution proportional 
to the elementary amplitudes p and P£° n as defined in Eq.(2). Their values, 

Pf p = 0.480 |q | GeV~ 2 and Pf n = 0.344 |q | GeV~ 2 have been taken from the 
P-wave low-energy theorems found in Bernard et al. (1996(1)) with q = /i (1 — 
/i) p — /i 2 m (1 — 5/i/4) k/2 (/i = M n /m and p is the nucleon three-momentum in 
the n°d COM system). E^ s has been evaluated using the Argonne V18, the Reid 
soft core (RSC),thc Nijmegen and the Paris potential, E s d s = 0.36 x 10 _3 /M 7r + 
with an uncertainty of 5E^ S = 0.05 x 10 _3 /M 7r + due to the various potentials 
used. The P-wave contribution amounts to a 3% correction to the one from 
the S-wave, i.e. it amounts to a minor correction. The sensitivity of the single- 
scattering contribution E^ s to the elementary neutron-7r° amplitude is given 

by 

E s d s = 0.36 - 0.38 • (2.13 - Ef + n ) , (13) 

all in units of lO^/M^. Consequently, for E]f + n = 0, one has E s d s = -0.45 
which is of opposite sign to the value based on the chiral perturbation theory 
prediction for Ey + n . If one were to use the empirical value for the proton ampli- 
tude, the single-scattering contribution would be somewhat reduced. 

• three body contribution 



12 Veronique Bernard 




Fig. 2. Graphs contributing at order q 4 to neutral pion photoproduction off deuteron. 
The circles denote an insertion from C^ N . 

To 0(q 3 ) only two diagrams contribute at threshold, see Beane et al. (1995). 
In one the photon couples to the pion in flight, the other is a seagull term 
which involves the charge exchange amplitude which is expected to dominate 
the single scattering contribution. To 0(q 4 ) few more graphs have to be added, 
these are shown in Fig. 2. These graphs do not involve any new unknown LEC. 
Furthermore there is no contribution from a possible four-fermion contact terms. 
One can therefore calculate E d in a parameter-free manner. One finds: 

E d = E s d s + Ef z + E d hA = 0.36 - 1.90 - 0.25 = -1.8 ± 0.2 (14) 

in units of 10~ 3 /M V +. The theoretical error is an educated guess. As expected 
E d b ' 3 is rather large. Note that tb, 4 is much smaller than tb, 3 which signals a 
good convergence. To see the sensitivity to the elementary neutron tt° amplitude, 
one sets the latter to zero and find Ed = —2.6 x 10~ 3 /M n + which is considerably 
different from the chiral theory prediction, Eq.(14). For other values of EqJ 1 , E d 
can be calculated from Eq.(13). Obviously, the sensitivity to the neutron am- 
plitude is sizeable and is not completely masked by the larger charge-exchange 
amplitude as it is often stated. On the experimental side, there exist one deter- 
mination of Ed, which is a reanalysis of older Saclay data by Argan et al. (1988) 
giving E d = -1.7 ± 0.2 x 10 _3 /M 7r+ in good agreement with the CHPT result. 
This number however has to be taken with care since the extraction of the em- 
pirical number relies on the input from the elementary proton amplitude to fix 
a normalization constant. A more precise experimental determination of E d has 
just began. Indeed recent data taking have been performed in Saskatoon with 
the preliminary result E d — —1.45 ± 0.09 x 10~ 3 /M„.+ , see these proceedings. 
This is somewhat smaller in magnitude than the CHPT value. The possible role 
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of higher order unitarity corrections has been discussed by Wilhelm at this work- 
shop. Also isospin breaking effects have to be included before one can draw any 
definite conclusion. Another source of information will come from an approved 
experiment at the Mainz Microton to measure the threshold cross section for 
coherent neutral pion electroproduction off deuterium at a photon virtuality of 
k 2 = — 0.075GeV 2 . It will certainly be necessary to have a CHPT calculation to 
compare with. 

6 Neutral pion electroproduction off protons 

Producing the pion with virtual photons offers further insight since one can ex- 
tract the longitudinal S-wave multipole L + and also novel P-wavc multipoles. 
The CHPT calculation proceeds in exactly the same way as for photoproduc- 
tion, see Bernard et al. (1996(3)). One has, however, some new counterterms. 
In addition to the expected form factor corrections, E(k 2 ) = L(k 2 ) oc M^k 2 br\ v 
(where 5r\ p is fixed from the knowledge of the isovector nuclcon radius) one 
has two new LECs, a 3 and a.4 for the S-wave multipoles which are such that by 
definition L-E~l + p(p= -k 2 /M 2 ), 



It turns out that these are strongly constrained by a soft-pion theorem, see 
Bernard et al. (1994). It implies that there is a strong correlation between the 
counterterms of E and L to order g 4 and furthermore that L ct is k 2 independent. 
To 0(q 4 ) with the P-dependence of L(k 2 ) coming solely from the Born and loop 
graphs one is unable to fit the existing data. The first corrections to the soft-pion 
constraint (03 + 04 = 0) away from the chiral limit have to be included. This 
induces terms of the type -Eq+> ~ 015 M 2 k 2 which are arising from terms in 

(5) 

the Lagrangian C^^ and are thus of higher order. These are the minimal terms 
one has to take to be able to describe the data at k 2 = —0.1 GeV 2 . Of course, 
there are other counterterms at this order. These, however, merely amount to 
quark mass renormalizations of the already considered fc 2 -independent counter 
terms and have therefore been set to zero. In Bernard et al. (1996(3)) a combined 
fit to the NIKEHF (e = 0.67), see van den Brink et al. (1995) and the MAMI 
(e = 0.582 and 0.885) , see Distler (1997), data at k 2 = -0.1 GeV 2 has been 
performed. Re Eo+, Re and the S-wave cross section ao = \Eo+\ 2 + cl\Lo+\ 2 
have been then predicted. As can be seen from Fig. 3 there is a very nice agree- 
ment between theory and experiment. One notices that Re E + has changed 
sign as compared to the photoproduction case, it shows the typical cups effect 
at the opening of the n + n threshold. In contrast, Re L 0+ is essentially energy- 
independent with a very small cusp. Note that ao is completely dominated by the 
L + multipole (dot-dashed line) since E + passes through zero at k 2 w —0.04 
GeV 2 . It should be however kept in mind that as one reaches values of k 2 of 
the order of —0.1 GeV 2 the one loop corrections are large so one should better 



E ct (k 2 ) 
L ct {k 2 ) 



eM 7I {(a 1 +a 2 ),M 2 -a 3 k 2 } , 

e {(ai + o 2 ) , M 2 -a 3 k 2 + a 4 (M 2 - k 2 



)} 



(15) 
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Fig. 3. Upper left pannel: Re Eo+, upper right pannel: Re Lo+ at k 2 = —0.1 GeV 2 . The 
diamonds are the NIKHEF data and the boxes the Mainz data. Lower part: S-wave 
cross section ao for e = 0.67 (solide line) in comparison to the data from Welch (boxes) 
and from van den Brink (cross). The dot-dashed line is the contribution from e|L 0+ | 2 . 

compare at lower virtualities. In Bernard et al. (1996(3)), many predictions for 
k 2 w -0.05 GeV 2 are given. At MAMI, data have been taken at k 2 = -0.05 
GeV 2 . The analysis is underway, see H. Merkel, these proceedings. Concerning 
the P-wave multipoles novel low-energy theorems have been derived to 0(q 3 ), 
sec Bernard et al. (1995). The combination M\ + — M\ — shows a rather good con- 
vergence. In the case of L\- and E\ + , an 0(q ) calculation is mandatory to 
have an idea of the next to leading order contributions. Note that these P-wave 
LETs have indeed been used in the analysis of the NIKHEF data, see van den 
Brink et al. (1995). 



7 Conclusions 

Since the MIT workshop considerable progress has been made: 

• charged pion photoproduction has been calculated to 0(q 4 ), 

• the energy dependence of E 0+ in the reaction -yp — > ir°p in the threshold 
region has been understood, 
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• a study of the sensitivity of ir° photoproduction off deuteron to the neutron 
amplitude has been performed showing a rather large one, 

• A fit to the data for n° electroproduction off protons at k 2 = —0.1 GeV 2 
has been performed allowing to make predictions at smaller virtualities, 

• another interesting process, which I had no time to discuss here, is double 
neutral pion photoproduction. It is found within an HBCHPT calculation to 
0(q i ) 1 see Bernard et al. (1996(4)), that this photoproduction channel is signifi- 
cantly enhanced close to threshold due to pion loops. The experimental analysis 
of the TAPS data is underway, see Stroher (1997). 

More precise low energy data as well as more refined calculations are still 
needed to further test the chiral dynamics of QCD, one of the most important 
question now being related to isospin violation in the pion-nucleon interaction 
and the role of electromagnetic corrections. Work along these lines has recently 
started, see Meifiner and Steininger (1997). 
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